Thermophysical properties for shock compressed polystyrene 
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We have performed quantum molecular dynamic simulations for warm dense polystyrene at high 
pressures. The principal Hugoniot up to 790 GPa is derived from wide range equation of states, 
where contributions from atomic ionizations are semiclassically determined. The optical conductiv- 
ity is calculated via the Kubo-Greenwood formula, from which the dc electrical conductivity and 
optical reflectivity are determined. The nonmetal-to-metal transition is identified by gradual decom- 
position of the polymer. Our results show good agreement with recent high precision laser-driven 
experiments. 

PACS numbers: 82.35.Lr, 51.30.-K, 52.27.Jt, 82.40.Fp, 71.15.Pd 



In high energy density (E/V^IO 11 J/m 3 ) physics 
(HEDP), pressure-induced response of materials, which 
can be probed through shock wave experiments jlH_|, is 
of crucial interest for inertial confinement fusion (ICF). 
Fundamental experiments in typical ICF designs, such 
as the investigation of hydrodynamic instabilities, have 
introduced plastic compounds as basic candidates to 
achieve high gain [5j . Furthermore, in recent fast ignitor 
experiments, the electrical properties of polymers have 
been often implicated for the transport of relativistic elec- 
tron beams in solid targets as compared to conducting 
materials [_| . Understanding the behavior of polymers at 
several megabar regime brings better insight into physi- 
cal processes in ICF. 

Some of the ablators in the indirect-drive mode of 
National Ignition Facility (NIF), will be made of glow- 
discharge polymer (GDP) with various levels of germa- 
nium doping (Ge-GDP) Due to the fact that no high 
pressure data exist for Ge-GDP, polystyrene (CH), which 
is closest in structure, is considered as a coarse indica- 
tor for shock timing simulations of NIF targets involv- 
ing such ablators [_|. The absolute equation of states 
(EOS) of polystyrene along the principal Hugoniot curve 
has been studied by using gas gun up to 50 GPa [__- 
[l2| and high energy laser-driven shock waves up to 4000 
GPa [13h15| . High energy dynamic compressions pro- 
duce the so-called warm dense matter (WDM), where si- 
multaneous dissociations, ionizations, and degenerations 
make the transition between condensed matter physics 
to plasma physics. The EOS and the relative properties, 
such as the Hugoniot curve, are important features in 
this context. Moreover, electrical properties such as the 
dielectric function are closely related to the dynamic con- 
ductivity, which determines the dc electrical conductiv- 
ity in the static limit. Then optical reflectance can also 
be extracted. All these quantities are indispensable in 
characterizing the unique behavior of shock compressed 
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polystyrene, especially for high pressures at, or exceed- 
ing, 1 Mbar. Recently, the development and applica- 
tion of quantum molecular dynamic (QMD) simulation 
16] techniques for WDM are available due to the enor- 
mous progress in computer capacity. QMD simulations, 
where quantum effects and correlations are systemati- 
cally treated, have conducted highly predictive results to 
describe WDM. 

In the present work, QMD simulations are applied to 
calculate a broad spectrum of thermophysical properties 
for shock compressed polystyrene. The self-consistent 
electronic structure calculation within density functional 
theory (DFT) yields the charge density distribution in 
the simulation supercell at every time step. The ion-ion 
pair correlation function (PCF), which is important for 
characterizing and identifying phase transitions, can be 
given by molecular dynamics run. The Hugoniot curve, 
which is derived from EOS data for a wide region of densi- 
ties and temperatures, is determined and compared with 
available experimental and theoretical results. As a start- 
ing point, we use the Kubo-Greenwood formula to eval- 
uate the dynamic conductivity a{uS), from which the dc 
conductivity, the dielectric function e{uS), and the optical 
reflectivity can be settled. 

We introduce ab initio plane wave code VASP [l__, [H| 
to perform QMD simulations. The elements of our 
calculations consist of a series of volume-fixed super- 
cells including N atoms, which are repeated periodically 
throughout the space. By involving Born-Oppcnheimcr 
approximation, electrons are quantum mechanically 
treated through plane-wave, finite-temperature DFT 
[l_j, where the electronic states are populated accord- 
ing to Fermi-Dirac distributions. Sufficient occupational 
bands are included in the overall calculations (the occu- 
pational number down to 10 -6 for electronic states are 
considered). The exchange-correlation functional is de- 
termined by generalized gradient approximation (GGA) 
with the parametrization of Perdew-Wang 91 [20]. The 
ion-electron interactions are represented by a projector 
augmented wave (PAW) pseudopotential [2l| . Isokinetic 
ensemble (NVT) is adopted in the present simulations, 
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FIG. 1: (Color online) The P-V (left panel) and T-P (right 
panel) Hugoniot curves of shocked polystyrene. Inset is the 
(u s , Up) diagram. Previous results are also shown for compar- 
ison. Experiments: gas-gun results (the gray filled squares); 
absolute measurements on NOVA (magenta diamonds); 

1311 high energy laser-driven measurements by Ozaki et al. 

3 EH (blue upward and downward triangles) and Barrios et 
al. (orange squares). [| Theories: QEOS and SESAME [H 
results are denoted by green dashed line and purple dotted 
line, respectively. 



where the ionic temperature Ti is controlled by Nose ther- 
mostat [22j . and the system is kept in local equilibrium 
by setting the electron (T e ) and ion (Tj) temperatures to 
be equal. 

The simulations have been done for 128 atoms (namely, 
eight CsHg units) in a supercell. A plane-wave cutoff en- 
ergy of 600.0 eV is selected, so that a convergence of bet- 
ter than 3% is secured. The Brillouin zone is sampled by 
r point and 3x3x8 Monkhorst-Pack scheme [23[ k points 
in molecular dynamic simulations and electronic struc- 
ture calculations, respectively, because EOS (conductiv- 
ity) can only be modified within 5% (15%) for the selec- 
tion of higher number of k points. The densities adopted 
in our simulations range from 1.05 g/cm 3 to 3.50 g/cm 3 
and temperatures between 300 and 50000 K, which high- 
light the regime of principal Hugoniot. All the dynamic 
simulations are lasted 6000 steps, and the time steps for 
the integrations of atomic motion are selected according 
to different densities (temperatures) [24j. Then, the sub- 
sequent 1000 step simulations are used to calculate EOS 
as running averages. 

The nature of high pressure behavior for warm dense 
polystyrene is dominated by a two-stage transition, that 
is, dissociations and ionizations. QMD simulations have 
been demonstrated to be powerful in describing as well 
as understanding chemical decompositions and electronic 
excitations of shock compressed materials, namely, the 
shocked atomic and electronic structures can be effec- 
tively treated within QMD. Importantly, however, we 
should address that the calculated EOS data from di- 
rect QMD simulations on shocked polystyrene should be 
corrected due to the following aspects: (i) At laser-driven 
temperatures as high as of 5~6 eV (~10 Mbar), the con- 



TABLE I: Polystyrene Hugoniot results from QMD simula- 
tions. 

p (g/cm 3 ) P (GPa) T (eV) u s (km/s) u p (km/s) 
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siderable effect of atomic ionizations on EOS should be 
appropriately included beyond QMD; (ii) QMD energy 
consists of kinetic and potential energy (for both ions 
and electrons), while, the well known zero-point vibra- 
tion energy (ZPVE), van der Waals energy (VDWE), and 
so as the demanded energy for phase transitions, which 
are particularly significant in WDM, are excluded; (iii) 
The pressure from QMD only describes the interaction 
contributions, and the ideal part for nonintcracting par- 
ticles are missed. To overcome these obstacles, in our 
present calculations for polystyrene, ZPVE, which is im- 
portant for determining the low-pressure internal energy 
for the reference state along the Hugoniot curve, is sim- 
ply added onto the QMD output. Condensed matter to 
plasma transition energy mainly comes from atomic ion- 
ization energy, which is also added onto QMD data, and 
the reference ionization degree is semiclassically deter- 
mined [25]. VDWE has been treated as small as neg- 
ligible. As for another important parameter — pressure, 
contributions from both interacting and noninteracting 
parts (ions and free electrons, separately) are included in 
the present EOS data. 

High precision EOS data are essential for understand- 
ing the electrical and optical properties of polystyrene 
under extreme conditions. The present EOS have been 
examined theoretically through the Rankine-Hugoniot 
(RH) equations, which follow from conservation of mass, 
momentum, and energy across the front of shock waves. 
The locus of points in (E 7 P, V)-sp&ce described by RH 
equations satisfy the following relations: 

E 1 -E = l{P 1 +P )(V -V ), (1) 
P 1 -P = p n u s u p , (2) 

V 1 = V (l-u p /us), (3) 

where subscripts and 1 present the initial and shocked 
state, and E, P, V denote internal energy, pressure, vol- 
ume, respectively. The u s and u p correspond to the shock 
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and mass velocities of the material behind the shock 
front. As the starting point along the Hugoniot, the ini- 
tial density is p =l-05 g/cm 3 , and the internal energy is 
Eq=— 87.44 kJ/g at a temperature of 300 K. Compared 
to the high pressure of shocked states along the Hugoniot, 
the initial pressure Pq can be treated approximately as 
zero. We use smooth functions to fit the internal energy 
and pressure in terms of temperature at sampled density, 
and derive Hugoniot point from Eq. [TJ The calculated 
Hugoniot data are listed in Table [TJ 

The principal Hugoniot curve is plotted in Fig. [TJ 
where previous theoretical and experimental results are 
also shown for comparison. The EOS of polystyrene has 
been previously probed by gas gun experiments, which 
have the advantage of high precision, but the pressure 
hardly exceeds 50 GPa. Recently, high energy laser- 
driven experiments have detected pressures up to 4000 
GPa. Meanwhile, however, large error bars were intro- 
duced at above 100 GPa [13, LLJ], and thus the use of 
low-precision EOS of the polymer to predict the behav- 
ior of NIF Ge-GDP ablator materials provides an un- 
acceptable uncertainty. Then, high precision EOS (up 
to 1000 GPa) for polystyrene was obtained by using a- 
quartz as an impedance- matching (IM) standard |8j|. As 
shown in Fig. [TJ previous data by Ozaki et al. fl4jl (IM 
with an aluminum standard) and Cauble et al. }13j (ab- 
solute data) show clearly stiffer behavior compared to 
those results with quartz standard [1, [lj| and SESAME 
model 26] . The possible reason for this stiff behavior is 
likely due to x-ray preheating of the samples, as has been 
stated by those authors. Then, thicker pushers and low- 
Z ablators were used in the newly measured data from 
Ozaki et al. [l5[ to reduce pre-heating of the samples. 
The experiments, where IM with a quartz standard were 
also used, show results that are closer to those by Bar- 
rios et al. In our QMD simulations, we find that direct 
calculated EOS without corrections mentioned above are 
only valid in the low pressure (P<100 GPa and T<1 eV) 
regime, beyond which no Hugoniot points can be found 
from the pure QMD data. For higher pressures, due to 
our introduction of corrections to direct QMD results, 
the calculated principal Hugoniot curve are greatly soft- 
ened, and the results show good agreements with those 
experiments where quartz standard was used. 

Temperature, as has been focused as one of the most 
important parameters in experiments, is difficult to be 
measured because of the uncertainty in determining the 
optical-intensity loss for ultraviolet part of the spectrum 
in adiabatic or isentropic shock compressions, especially 
for the temperature exceeding several eV [27| ■ QMD sim- 
ulations can provide efficient predictions for shock tem- 
peratures. As has been shown in Fig. [JJ (right panel), our 
calculated Hugoniot temperatures are accordant with ex- 
perimental results up to 2 eV, and discrepancy emerges 
at higher temperatures. Beyond 2 eV, the predicted tem- 
peratures from both experiments and SESAME model at 
a given pressure are higher compared to our calculations. 

We also examine the structural transition of 
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FIG. 2: (Color online) Pair correlation functions for four dif- 
ferent points along the principal Hugoniot curve, (a) ~ (d) 
correspond to the starting point, 53.48 GPa, 74.54GPa, and 
128.56 GPa, respectively. C-H, C-C, and H-H bonds are de- 
noted by blue, black, and red lines, respectively. Insets are 
the sampled atomic structure and charge density distribution 
at each state (gray balls for carbon and white balls for hydro- 
gen). 



polystyrene under shock compressions by using PCF, 
which is evaluated at equilibrium during the molecular 
dynamics simulations. Along the Hugoniot, four points 
of PCF are shown in Fig. [2] At the initial state [see Fig. 
[2a)], the ideal condensed polymer phase is indicated by 
one peak for C-H bond, meanwhile, two peaks (7r-type 
and sp 3 -type) for C-C bond. With the increase of pres- 
sure (~20 GPa), phenyl decomposes, which is indicated 
by that the two C-C peaks merge together (a transverse 
from 7r bond to sp 2 , sp 3 like bond). At this stage dia- 
mondlikc carbon nanoparticles (with defects) are formed. 
The present phenomenon agrees well with the behavior 
of shock compressed benzene [HI, [29} . Further increase of 
pressure produces molecular and atomic hydrogen, which 
is indicated by the peak in H-H PCF around 0.75 A [Fig. 
[2jb)]. Decomposition of hydrocarbons continues until 
128 GPa, where PCF shows no obvious evidence for the 
existence of C-H bond. Whereas, for higher pressures, 
the PCFs do not show visible difference, and are thus 
not presented here. SESAME model suggested that the 
softening behavior of Hugoniot is caused by the chem- 
ical decomposition of the polymer at 200 ~ 400 GPa. 
Whereas, our calculations show that the decomposition 
pressure lies between 20 and 128 GPa, which is much 
lower than that of SESAME. Note that the chemical de- 
compositions do not contribute much to the softening 
behavior, but the ionizations really do. 

Having clarified the EOS, let us turn now to study the 
electrical and optical properties for shocked polystyrene. 
Based on Kubo-Greenwood formula, the electrical and 
optical properties can be obtained as described in Ref. 
28]. In order to get converged results, 30 indepen- 
dent snapshots, which are selected during one molecu- 
lar dynamics simulation at given conditions, are picked 
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FIG. 3: (Color online) Dc conductivity (left panel) and op- 
tical reflectivity at the wavelength of 532 nm (right panel). 
Previous measurements by Koenig et al. (green squares) [30] 
and Barrios et al. (upward triangles) [8J are also plotted. 

up to calculate dynamic conductivity as running aver- 
ages. The dc conductivity ddc, which follows from the 
static limit u>— >0 of <ti(lo) is then evaluated and plot- 
ted in Fig. [3] (left panel) as a function of shock ve- 
locity. Initially, er^c increases rapidly with shock ve- 
locity up to 14 km/s (128 GPa) towards the formation 
of metallic state of polystyrene. Then, one can find 
from Fig. [3] that ddc keeps almost invariant and the 
shock compressed polystyrene maintains its metallic be- 
havior. The onset of metallization (<7,j! c >10 3 f2 — 1 cm — •*•) is 
observed at around 10 km/s (—50 GPa), where dissocia- 
tion of the polymer and formations of hydrogen (molecu- 
lar and atomic species) govern the characteristics of warm 
dense polystyrene. We stress here that the nonmetal to 
metal transition is induced by gradual chemical decom- 
positions and thermal activations of the electronic state, 
instead of atomic ionization, which is not observed un- 
til 128 GPa with respect to the charge density distri- 
bution in the QMD simulations. In the same pressure 
range (20—200 GPa), we observe larger Udc compared to 
experimental measurement (30| . This overestimation of 
Udc values could be attributed to the use of DFT-based 
molecular dynamics, which is known to underestimate 
band gaps in many systems. 

Optical reflectance, from which emissivity can be de- 
rived, is of great interest in experiment for determin- 



ing shock temperature. Along the Hugoniot, we show 
in Fig. [3] (right panel) the variation of optical reflectance 
at 532 nm as a function of shock velocity. Optical prop- 
erties of polystyrene have been experimentally studied 
by Koenig et al. [H with u s of 11-16 km/s (80-170 
GPa), and steadily increasing reflectivities reaching val- 
ues up to —50 % were observed. The results by Ozaki et 
al. indicated reflectivity from 16% to 42% in u s range of 
22-27 km/s (300-500 GPa) [H|. Recent experiments by 
Barrios et al. [|| suggested a drastic increase in reflectiv- 
ity around the Hugoniot pressure of 100 GPa, followed 
by saturated value of 40% around 250—300 GPa. Dis- 
crepancies in reflectivities among experiments have been 
claimed to be arisen from probe-beam stability or from 
differences in diagnostic configurations. Along the Hugo- 
niot, our QMD calculations provide the general feature 
for the optical reflectance — steep increase (from 15% to 
60 %) followed by saturation (at —128 GPa). The change 
in optical reflectivity with pressure can be interpreted as 
the gradual insulator-conductor transition at above 20 
GPa, where the atoms strongly fluctuate with neighbors 
to dense, partially ionized plasma at high pressures above 
128 GPa. 

In summary, we have carried out QMD simulations 
to study the thermophysical properties for warm dense 
polystyrene. The Hugoniot EOS data of polystyrene up 
to —790 GPa, which is in agreement with dynamic ex- 
periments in a wide range of shock conditions, has been 
evaluated through QMD calculations and corrected by 
taking into account the atomic ionization. A two-stage 
(dissociation and ionization) transition governs the char- 
acteristic of polystyrene under extreme conditions. Con- 
tribution from chemical decomposition demonstrates the 
steep increase of Ode and optical reflectance observed at 
20—128 GPa. While, soften behavior of the Hugoniot is 
dominated by atomic ionizations for higher pressure. 
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